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Possible explanations of the recently reported anomalous increase of the eccentricity of the 
lunar orbit are sought in terms of classical Newtonian mechanics, general relativity, and long- 
range modifications of gravity. 

Anderson and Nieto, in a recent revie of some astrometric anomalies detected in the 
solar system by some independent groups, mentioned also an anomalous secular increase of the 
eccentricity e of the orbit of the Moon 

Ws = (9 ± 3) X 10-12 yr-i (1) 

based on an analysis of a long LLR data record spanning 38.7 yr performed by Williams and 
BoggiP with the dynamical force models of the DE421 ephemeride^'^ including all the known 
relevant Newtonian and Einsteinian effects. Notice that Eq. [1] is statistically significant at a 
3a"— level. The first accoun^i^ of this effect appeared in 2001 by Williams et al, who gave an 
extensive discussion of the state-of-the-art in modeling the tidal dissipation in both the Earth and 
the Moon. Later, Williams and Dicke;^, relying upon the 2001 stud;^, released an anomalous 
eccentricity rate as large as emeas = (1-6 ± 0.5) x 10-^^ y^-^- Anderson and Niet(P commented 
that Eq. [T]is not compatible with present, standard knowledge of the dissipative processes in the 
interiors of both the Earth and Moon, which were, actually, modeled by Williams and BoggP^. 

Naive, dimensional evaluations of the effect caused on e by an additional anomalous accel- 
eration A can be made by noticing that 

ec^—, (2) 
na 

with 

na = 1.0 X 10^ m s^^ = 3.2 x 10^° m yr^^ (3) 

for the geocentric orbit of the Moon, whose mass is denoted as m. In it, a is the orbital 
semimajor axis, while n = ^yjt/a^ is the Keplerian mean motion in which // = GM{1 + m/M) 
is the gravitational parameter of the Earth-Moon system: G is the Newtonian constant of 
gravitation and M is the mass of the Earth. It turns out that an extra-acceleration as large as 

A ~ 3 X 10^1^ m s-2 = 0.3 m yr'^ (4) 



would satisfy Eq. [H In fact, a mere order-of-magnitude analysis based on Eq. [2] would be 
inadequate to infer meaningful conclusions: finding simply that this or that dynamical effect 



induces an extra-acceleration of the right order of magnitude may be highly misleading. Indeed, 
exact calculations of the secular variation of e caused by such putative promising candidate 
extra-accelerations A must be performed with standard perturbative techniques in order to 
check if they, actually, cause an averaged non-zero change of the eccentricity. Moreover, it may 
well happen, in principle, that the resulting analytical expression for (e) retains multiplicative 
factors l/e^,j = 1,2,3,... or e^,j = 1,2,3... which would notably alter the size of the found 
non-zero secular change of the eccentricity with respect to the expected values according to 
Eq.El 

It is well known that a variety of theoretical paradigms^ allow for Yukawa-like deviation^ 
from the usual Newtonian inverse-square law of gravitation. The Yukawa-type correction to the 
Newtonian gravitational potential ?7n = ~/^/'') where /i = GM is the gravitational parameter 
of the central body which acts as source of the supposedly modified gravitational field, is 

C/v = -^exp(-n, (5) 



X 

in which ^oo is the gravitational parameter evaluated at distances r much larger than the scale 
length X. In order to compute the long-term effects of Eg. [5lon the eccentricity of a test particle 
it is convenient to adopt the Lagrange perturbative scheme?^. In such a framework, the equation 
for the long-term variation of e i m 

/de\_ 1 / l-e^ W 1 dTZ dn\ 

where uj is the argument of pericenter, A4 is the mean anomaly of the test particle, and TZ 
denotes the average of the perturbing potential over one orbital revolution. In the case of a 
Yukawa-type perturbation, Eq. [5] yields 



g^Uoo exp {-j) ^ /ae\ 



= - - [^) , (7) 

where Iq{x) is the modified Bessel function of the first kind Iq{x) for q = 0. An inspection 
of Eq. [6] and Eq. [7| immediately tells us that there is no secular variation of e caused by an 
anomalous Yukawa-type perturbation. 

The size of the general relativistic Lense-Thirrin^I^ acceleration experienced by the Moon 
because of the Earth's angular momentunff2! g = 5_8g x 10^'^ kg m^ s^^ is just 

2vGS 



2n3 



1.6 X 10"^^ m s~2 = 0.16 m yr"^ (8) 



i.e. close to Eq. HI On the other hand, it is well known that the Lense-Thirring effect does not 
cause long-term variations of the eccentricity. Indeed, the integrated shift of e from an initial 
epoch corresponding to /o to a generic time corresponding to / iff^ 

2GS cos / (cos f — cos fn) ,„s 
Ae = ^ ■' — (9) 

c^na^v 1 — 

in which /' is the inclination of the Moon's orbit with respect to the Earth's equator and / is 
the true anomaly. From Eq. [9] it straightforwardly follows that after one orbital revolution, i.e. 
for / — 7- /o + 27r, the long-term gravitomagnetic shift of e vanishes. 

A promising candidate for explaining the anomalous increase of the lunar eccentricity is, 
at least in principle, a trans-Plutonian massive body X of planetary size located in the remote 
peripheries of the solar system. Indeed, the perturbation induced by it would, actually, cause a 



non-vanishing long-term variation of e. Moreover, since it depends on the spatial position of X 
in the sky and on its tidal parameter 



= (10) 

where mx and dx are the mass and the distance of X, respectively, it may happen that a suitable 
combination of them is able to reproduce Eq. [TJ Let us recall that, in general, the perturbing 
potential felt by a test particle orbiting a central body due to a very distant, pointlike mass can 
be cast into the following quadrupolar form 



Ux- — 



-3 ff • 



(11) 



where I = {Ix, ly,lz} is a unit vector directed towards X determining its position in the sky. In 
Eq. [TT] r = {x, y, z} is the geocentric position vector of the perturbed particle, which, in the 
present case, is the Moon. loricP^has recently shown that the average of Eq. [TTjover one orbital 
revolution of the particle is 

{Ux) = '^u[e,I,n,oj-j), (12) 

where U (^e, I,Q,uj;tj is a complicated function of its argumentiP: is the longitude of the 
ascending node and / is the inclination of the lunar orbit to the ecliptic. In the integration 
yielding Eq. [T2]'I was kept fixed over one orbital revolution of the Moon, as it is reasonable given 
the assumed large distance of X with respect to it. Eq. [6l applied to Eq. [I2l straightforwardly 
yields 

{e) = '-^^^^^f^£(l,n,u;i). (13) 



16n 

Also £ [I,Q,,uj;l\ is an involved function of the orientation of the lunar orbit in space and of the 



position of X in the sky^. Actually, the expectations concerning X are doomed to fade away. 
Indeed, apart from the modulation introduced by the presence of the time- varying I,u} and O 
in Eq. [T3l the values for the tidal parameter which would allow to obtain Eq. [T] are too large for 
all the conceivable positions {/3x, Ax} of X in the sky. This can easily be checked by keeping u 
and fixed at their J2000.0 values as a first approximation. Indeed, lorio^^ showed th at the 
physical and orbital features of X postulated by two recent plausible theoretical scenarioP^I^ 
for X would induce long-term variations of the lunar eccentricity much smaller than Eq. [TJ 
Conversely, it turns out that a tidal parameter as large as 

/Cx = 4.46 X 10-2^ s-2 (14) 

would yield the result of Eq. [TJ Actually, Eq. [HJ is totally unacceptable since it corresponds to 
distances of X as absurdly small as dx = 30 au for a terrestrial body, and dx = 200 au for a 
Jovian mass. 

An empirical explanation of Eq. [TJcan be found by assuming that, in addition to the usual 
Newtonian inverse-square law for the gravitational acceleration imparted to a test particle by a 
central body orbited by it, there is also a small radial extra-acceleration of the form 

A = kHoVr. (15) 

In it is a positive numerical parameter of the order of unity to be determined from the 
observations, Hq = (73.8 ± 2.4) km s'^ Mpc'^ = (7.47 ± 0.24) x lO^^^ yi'^ is the Hubble 



parameter at the present epocllJ^, defined in terms of the time- varying cosmological scahng factor 
S{t) as Hq = S/S , and Vr is the component of the velocity vector v of the test particle's proper 
motion about the central body along the common radial direction. Indeed, a straightforward 
application of the Gauss perturbative equation for e to Eq. [15] yields 



e2 



(e) = kHo ^— . (16) 

Since eMoon = 0.0647, Eq. [16] can reproduce Eq. [1] for 2.5 ^ k < 5. Here we do not intend to 
speculate too much about possible viable physical mechanisms yielding the extra-acceleration 
of Eq. [15] It might be argued that, reasoning within a cosmological framework, the Hubble 
law may give Eq. [15] for A: = 1 if the proper motion of the particle about the central mass is 
taken into account in addition to its purely cosmological recession which, instead, yields the 
well-known local extra-acceleration of tidal type j4cosmoi = —QoHqV, where qo is the deceleration 
parameter at the present epoch. 
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